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OfT-diagonal geometric phases have been developed in order to provide information of the geometry 
of paths that connect noninterfering quantal states. We propose a kinematic approach to off-diagonal 
geometric phases for pure and mixed states. We further extend the mixed state concept proposed in 
[Phys. Rev. Lett. 90, 050403 (2003)] to degenerate density operators. The first and second order 
off-diagonal geometric phases are analyzed for unitarily evolving pairs of pseudopure states. 
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I. INTRODUCTION 

Pancharatnam's p]| work on the geometric phase in in- 
terference of light waves in distinct states of polarization 
plays a prominent role among early anticipations of the 
concept of geometric phase. It found its quantal coun- 
terpart when Berry 2] discovered geometric phase fac- 
tors accompanying cyclic adiabatic changes. Since then 
there has been an immense interest in holonomy effects in 
quantum mechanics, which has led to many generaliza- 
tions of the notion of geometric phase. The extension to 
nonadiabatic cyclic evolution was developed by Aharonov 
and Anandan 3, 4| . Samuel and Bhandari Q generalized 
the pure state geometric phase further by extending it to 
noncyclic evolution and sequential projection measure- 
ments. Mukunda and Simon put forward a kinematic 
approach to the geometric phase. Further generalizations 
and refinements, by relaxing the conditions of adiabatic- 
ity, unitarity, and the cyclic nature of the evolution, have 
since been carried out [3, |3| . 

The concept of geometric phases for pure states has 
been extended to mixed states. This was first addressed 
by Uhlmann 0, [l^ within the mathematical context of 
purification. Sjoqvist et al. introduced an alternative 
definition of geometric phase for density operators with 
all the nonzero eigenvalues being nondegenerate, based 
upon the experimental context of quantum interferome- 
try. Singh et al. |l2i | gave a kinematic description of the 
mixed state geometric phase in Ref. [Tll | and extended 
it to degenerate density operators. The generalization 
to nonunitary evolution has also been addressed in Refs. 
p3.. .14 jjl and experimental test of the phase concept 
in Ref. has been carried out using nuclear magnetic 



resonance [T^ . 

All the above notions of geometric phase break down in 
cases where the interference visibility between the initial 
and final states vanishes leading to an interesting nodal 
point structure in the experimental parameter space 
that could be monitored in a history-dependent manner 
'it', "ll 111 IS mill- ThisprobW may be overcome by 
introducing the concept of off-diagonal geometric phases, 
as was first put forward for pure states in adiabatic evo- 
lution by Manini and Pistolesi '2^ . Mukunda et al. pl | 
removed the adiabaticity assumption of Ref. [23j and 
Hasegawa et al. |2^ |2a| verified the off-diagonal pure 
state geometric phase by a neutron experiment. The con- 
cept of off-diagonal geom etric phase has been extended 
to mixed states |23, represented by nondegenerate 
density operators, as well as to the Uhlmann holonomy 
|2^ . These off-diagonal mixed state geometric phases 
contain information about the geometry of state space 
along the path connecting pairs of density operators, 
when the standard mixed state geometric phases in Refs. 
H lift ITlj are undefined. In this paper, we propose a 
kinematic approach to the off-diagonal mixed state geo- 
metric phase in Refs. [27l f28l | and further extend it to 
the degenerate case. 



II. KINEMATIC APPROACH TO 
OFF-DIAGONAL GEOMETRIC PHASES 

Consider an N dimensional quantum system s, the 
pure states of which represented by vectors of Hilbert 
space Tis. Let . . . jIiPn) denote an arbitrary or- 
thonormal basis of and {UHt)\t G [0, r], J/H (0) = /} 
be a one-parameter family of unitarities fulfilling the par- 
allel transport conditions 



{ijk\U\\^{t)uHt)\ijk) =0, fc = 1,...,7V 



(1) 
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relative to the chosen basis. For any pair of basis vec- 
tors and I'f/'ja) such that {ipji\U^^{T)\^pj^) ^ and 



2 



(V'j2|t^"WIV'ji) + 0, the quantity 



(2) 



(V',jC/ll(r)|^,,)(V',jC/"(T)|V,J 



(2) 



where $[z] = z/|z| for any nonzero complex num- 
ber z, is gauge invariant and therefore a property 



of the pair of paths Cji,C 



32 



[0/ 



IV'ii WMV'js W)('0i2(i)l in projective Hilbert 

space. The set {7) J,} constitutes the second order off- 
diagonal pure state geometric phase factors "23]. It may 
be extended by considering / < N states yielding the Ith 
order off-diagonal pure state geometric phase factors 



{I) 



n'J'[(V',Jt/ll(r)|^,„,J 



(3) 



where IV'ji+i) — IV'ji) and (V'ja|f^"(''')|V'ja+i) for all 
a = 1, . . .1. 

Notice that there is an equivalence set S of unitarities 
U(t) that all realize {Ck}, namely those of the form 



N 



(4) 



k=l 



where U{t) G S and 0k{t) are realtime-dependent param- 
eters such that dk{0) = 0. We may in particular identify 
U^t) e 5 by substituting U^t) = U(t) into Eq. id, so 
that we obtain 



{t)^i {i^k\U\t')U{t')\^k)dt' (5) 



and 



^ rt t , ■ 

[/II (t) = {/(t)^e-/o<^'=l^'(* mm (6) 



k=l 



Inserting this expression into Eq. Q for [/"(t), we ob- 
tain a kinematic expression for the I = 2 off-diagonal 
geometric phase factors as 



(2) 



- $[(^,jC/(r)|V,.)(V',Jt/(r)|V,J 



(7) 



In the above expression, the unitary operator U{t) G S 

need not satisfy the parallel transport conditions. One 

(2) (2) 
may verify that Tj,-jJ^(^^^ = 7]ijjc/(t) for any choice of 

dkit). 

An interesting aspect of Eq. ([TJl is that in the qubit 
case with 1-01) |0) and 1-02) = it follows that 
(0|C/t(t)[/(t)|0) = -(l|C/t(i)f/(t)|l) for any SU(2) oper- 
ation, leading to 



(2) 

7oi 



(0|C/(r)|l)(l|L/(r)|0) 
(0|C/ll(r)|l)(l|t/ll(r)|0) 



-1, 



(8) 



whenever (0|?7(t)|1) 7^ and (l|[/(r)|0) ^ 0. Thus, the 
TT shift in the qubit case is independent of the fulfillment 
of the parallel transport condition Eq. (Q; a result that 
has been experimentally verified for neutron spin |25ll26| . 

Similarly, for the higher order off-diagonal geometric 
phases ,2^, we obtain 



a=l 

xe 



n<I'[(^,Jt/(r)|V',„+J 



upon substitution of Eq. © into Eq. Once again, 

one may verify that 7S..jj£?(t) = lj[\..j,\u{t)- 

We now put forward a kinematic approach to the off- 
diagonal geometric phases of nondegenerate mixed states 
proposed in Refs. JL^il- Let 



pi = Xi\-ipi) {ipi 



AAr|0Ar)(V'Ar| 



(10) 



where for all nonzero eigenvalues we have ^ Xi^^k- 
Furthermore, introduce the unitarity 

W = \yJl){4'N\ + \^n){^N~i\ + ■■■ + |^2)(0l| (11) 



in terms of which any pair of members pj^, pj^, ji 7^ j2, 
of the set 



Pn^w^-'p,(w^y 



= l,...,iV (12) 



are unitarily connected and do not interfere since 
Tr(VF^^^-'^Pjj) vanishes. In analogy with the pure state 
case, the density operators pn constitute a set of states 
in terms of which we may define the off-diagonal mixed 
state geometric phase factors as [23, l2^ 



^(0 



Tr 



(13) 



as long as the trace in the argument of <I> does not vanish. 
Here, ?7"(i) satisfies the parallel transport conditions in 
Eq. and I < N. Apparently 7^^^^ are the standard 
mixed state geometric phase factors associated with the 
unitary paths A4n t G [0,t] Pn{t) in state space. 
Furthermore, it has been demonstrated [23, a phys- 
ical scenario for the / ~ 2 case in terms of two-particle 
interferometry. 

As the nonzero eigenvalues of all mixed states under 
consideration are nondegenerate, the parallel transport 
conditions are still given by Eq. CJ and U^t) by Eq. 
®. Substituting Eq. ® into Eq. (^21 7 we obtain the 
kinematic expression for the off-diagonal geometric phase 
factors for mixed states under evolution U (t) as 
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0,(0 = $ 

'Pii ■■■Pii 



N 



(V-., I C/(t) I ) (V-., I f/(T) I V-'.s > • • ■ (V-.™ I C/(t) I V'^, > 



ii ,...,ii=l 



X exp 



a=l 



(14) 



where A_ 



A 



N-p 



, p = 0, ...,7V — 1. One may verify t Pn{t) with p„(0) = p„, could be taken as 



that the phase factors ^pj^...pj^ are gauge invariant in that 
they are independent of the choice of U{t) e S. 



= 0, A: = 1, 



.K. (17) 



III. OFF-DIAGONAL GEOMETRIC PHASES 
FOR DEGENERATE DENSITY OPERATORS 

We now generahze the concept of off-diagonal geomet- 
ric phases of nondegenerate mixed states to degenerate 
mixed states. Let 



Pi — M^ia 



(15) 



where the nonzero eigenvalue Afc, k = 1, . . . , K < N , is 
TOfe fold degenerate and P]^''^ the corresponding projec- 
tor of rank mfe. The choice \ipi), . . . , \iPn) of orthonormal 
Hilbert space basis and concomitant permutation unitar- 

ity W = \^i){iPn\ + |?/'Jv)(V'w-i| + • • • + |'/'2)(V'i|, defines 
a set of nonintcrfcring density operators 



\ p(™i) 



\ p(™/<) „ 



where P^!^"' = ^^"-ipg'^VKt)"-!, n 1, , 



,N, (16) 



.N. Note 



thatPS-^^^P-r forn^n'. 

The parallel transport conditions for a unitary path 



In terms of the orthonormal Hilbert space basis 
[■i/ji), . . . , \i^N), Eq. p7|l is equivalent to the more familiar 
parallel transport conditions |3| {"^tiW}} (t)(j\{t)\4'v) = 
0, V pairs € Pjvk^^ ^^'^ reduces to those of Ref. 

[lll | in the nondegenerate case where K = N. 

Now, suppose there is a one-parameter family of uni- 
tarities {U{t)\t e [0,t],C/(0) = /} defining the paths 

Mn ■■ t e [0,t] ^ p„(<) = C/(i)p„C/t(t), 71 = l,...,iV. 

Then, for each p„ , there is an equivalence set 5„ of uni- 
tarities Unit) that all realize namely those of the 

form 



Unit) = Uit)Vnit) 



(18) 



with 



K(t) = a„;i(i) + ... + a«;i^(0, (19) 
where a„;fe(i) is unitary on the fcth degenerate sub- 
space, i.e., a„;fe(i)Q;^(t) = a,\.fe(i)a„;fe(i) = Pj^^^K and 
an;fe(0) = P!i"k''^ ■ We can identify C/,'|(t) G 5„ by substi- 
tuting Ujiit) = C/„(t) into Eq. lO, and obtain 



a. 



^ir'Texp 



P 



u\t')uit')Pll''Ut' 



Pi 



(mfc) 



(20) 



where T denotes time ordering. The parallel trans- 
porting unitary operators for p„ may be expressed 
as C/|(i) = Uit)Vnit) with supplementary operators 
vHt)^al,(t) + ... + al^{t). 

In the presence of degenerate subspaces, some of the 
parallel transporting unitarities u}i{t) must be differ- 
ent from each other. Thus, unlike the nondegenerate 
case, where a common parallel transporting J7"(t) can 
always be found, such a common operator parallel trans- 
porting all the states pn does not exist in the degen- 



erate case; we need to take f/j' (i), . . . , u\it) to parallel 



transport pi, 

-II 



, Pat, respectively. Substituting U^j^it) 



Uit)Vlit), a = 1,. into Eq. (O for the ath C/ll(i), 
we obtain the off-diagonal geometric phase factors of the 
degenerate mixed states as 



^(0 

'Pii ■■■Pii 



TA\{U{T)vlir) 



\a=l 



(21) 



for nonvanishing trace in the argument of $. The phase 
value calculated in Eq. (|21|) is manifestly gauge invariant 
under choice of the set {C/„(f)}, which shows that it is 
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only dependent upon the paths traced by the states. We 
may verify this point by repeating the demonstration of 
Ref. ^12]. 

Again, once the paths of states are defined by some uni- 
tary operator U{t), the supplementary operators vj[{t), 

a — 1,.. .J, and therefore C/j'^(t), may be different for 
each pj^ . However, if all the states pj^ have the same 
degeneracy structure [sij . they share a common unitary 
operator U^t). For a general set of states pj^, . . . ,pj^, 
whether a common unitary parallel transporting opera- 
tor exists depends both upon the states and the unitarity 
U{t). In the case where a common parallel transport- 
ing unitary operator does not exist, two or more unitary 
operators are necessary to calculate the off-diagonal ge- 
ometric phases. For this reason, experimental tests of 
off-diagonal geometric phases may be difficult to imple- 
ment in such cases. 



IV. EXAMPLE: PSEUDOPURE STATES 

Let us illustrate the above for pseudopure states, which 
take the form 

Pn = ^-^I + e\n)(n\, n=l,...,N>2 (22) 

with < e < 1, and \n) any p— qubit state. Such states 
have a single nondegenerate eigenvector \n) with eigen- 
value (l + {N — l)e)/iV and an — 1 fold degenerate 
eigenvalue ( 1 — e) /N with eigenprojector I—\n){n\. Pseu- 
dopure states appear generically as input states to stan- 
dard liquid-state NMR quantum computers, in case of 
which N = p being the number of nuclear spin qubits 



Now, we wish to compute the two lowest order {1 = 1 
and / = 2) geometric phases for the pseudopure states pn 
and pm with (n|m) = 0. Let 

U{t)=I-Pn^ + Unm{t), (23) 

where Pnm is the projection operator onto the sub- 
space spanned by {\n),\m)} and Un„iit)Ul^{t) = 
Ulmi^Wnmit) = Pnm- In the casc of liquid-state NMR 
quantum protocols, one may for example take 

\n) = |O)i0...(g)|O)g(8)...(^|O)p, 
|to) = |O)i0...(8)|I)5(g)...(^|O)p, (24) 

where 1 < q < p so that U{t) represents single-qubit 
operations on the qth nuclear spin. Using Eq. (|19(l . we 
obtain the parallel transporting unitarities 

C^l = Ul=I- P„™ + Un,n\n){n\ 

xe-/o<"l^"'"(*')^-(*')l">^*'+C/_|m)(m| 

^g-/J(m|C/+„(t')C/„™(t')|m>dt' (^25) 

Thus, due to the specific form of U, there is a single 
unitarity parallel transporting both /?„ and This 
property may be advantageous in experimental tests of 
the I = 1 and / = 2 geometric phases for p„ and pm- Note 
that a separate treatment of the random mixture (e = 0) 

yields Uji = Ul = I. 

Let us first analyze the / = 1 geometric phase factors 
for Pn and pm in the e > case. Inserting Eqs. (j22|l and 
(UHl into Eq. (EU yields 



(7W)* ^ $ J(iv _ 2)(1 - e) + 7j{l + {N- l)e)e-*^V2 + ^(1 _ 
exp — z arctan 



TjNe sin § 



{N - 2)(1 - e) + 77(2 + (A^ - 2)e) cos § 



(26) 



whenever the phase factor is defined. Here, fl is the 
geodesically closed solid angle on the Bloch sphere de- 
fined by projection of the two-dimensional subspace 
spanned by the pure state components \n) and |m) and 
rj = \{n\Unm\n)\ = |(m'|t',irn|m')| is the pure state visi- 
bility. In the absence of noise (e = 1) the I = 1 geo- 
metric phases become undefined iff i] = 0, corresponding 
to the case where |n) and |m) are interchanged. When 
77 ^ 0, we obtain the standard pure state geometric 
phases j^,]^ = (7p™)* — e~'^/^. In the general case 
(0 < e < 1), the nodal points, defined as those points 
in parameter space (77, e, Q) where j^]^ and "fp]^ are un- 



defined, are obtained as solutions of 

(^(A^ - 2) (1 - e) + 7,(2 + (A - 2)e) cos 

-1-772 A^^e^ sin^ ^ 0- (27) 

For N = 2 (nondegenerate case), the nodal points are 
those where 77 = 0, which is consistent with Refs. [2^12^ . 

and when 77 7^ we have 

^i'J = UllT = exp f-t arctan (etan \ , (28) 
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which is consistent with Refs. 

[Hill by identifying e as 
the length of the Bloch vector. For TV > 3 a necessary 
condition for a nodal point is Vl/2 = (2fc + l)7r, k integer, 
and 77 ^ 0. For these solid angles, Eq. (|27|l reduces to 



7? 



{N-2){l-e) 
2 + {N ~2)e ■ 



(29) 



Physical solutions (i.e., fulfilling ?7 < 1) of this equation 
exist provided 



e > 



2(iV- 2) 



(30) 



Thus, for N >b, there is a lower bound on the amount of 
noise for the existence of nodal points, i.e., in too noisy 
input states the I — 1 geometric phases are always well- 
defined for this U . 

Let us now turn to the I — 2 geometric phase for p„ and 



(2) 



(2) 

7p™p„ 



which is due to the 



Pm ■ First notice that 7p' 
symmetry of trace under cyclic permutations. Explicit 
calculation yields 



^(2) 



[N - 2)(1 -e) + {2 + {N~ 2)e)(-l + jf) 
+27f [l - e){l + {N - l)e) cosd. (31) 



Thus, TpJpm = ill changing sign across a nodal surface 
in the parameter space (e,ry, fi). This nodal surface is 
defined by vanishing argument of $ in Eq. H31|l . i.e., the 
solution of 

2 _ 2(A^ - 2)e - iV + 4 



2 + (iV - 2)e + 2y^(l-e)(l + (Ar-l)e) cos 

(32) 

If cos 57 > 0, solutions fulfilling r/ < 1 always exist for 
all N as long as Eq. (|30|l is satisfied. In the case where 
cosil < 0, however, < ry^ < 1 yields 



iV-4 
2(iV- 2; 



< e < 



(7V-2)2 -4cos2r2 
4(iV - 1) cos2 rj + (iV - 2)2 



(33) 



Thus, there are nodal points only for the amount of noise 
satisfying the above equation. In the limiting case N — 2, 
the upper bound on e is negative and thus no nodal points 
exist, which is consistent with Refs. 27, 28] . 

(21 

Now we wish to check 7p„p„, at the nodal points of 

hpll)*- For = 2 we obtain that the nodal 
points for I — 1 and I = 2 never coincide. For N > 3, 
insert ri/2 — {2k + l)7r, k integer, into Eq. H32f) and 
eliminate e by using Eq. (|29|l . yielding 



7^^ 



77-1 



2772 



iV- 2 



v/r;(7V - 2 - 77) = /(7;, N) - 0, (34) 



the solutions of which are the common nodal points of 



(2) 

7p„P„ 



and ^p]} ■ In Fig. ^we have depicted f{ri,N) for 



N = 3,4,5,6 and we conclude that there are coinciding 
nodal points for N > 2. For these rare parameter values 
one has to resort to Z > 3 geometric phases to retain 
some information about the geometry of the path in state 
space. 




FIG. 1: Common nodal points of 7p^'p„ and jp^J 



* for 



pseudopure states pn and pm are given by vanishing /{rj, N), 
which is depicted for A'^ = 3, 4, 5, 6. Thus, except for the qubit 
case {N = 2), there may arise situations where the 1 = 1 and 
I = 2 geometric phases are undefined simultaneously. 

(2) 

To measure 7pip2 we add an ancillary system a and 
consider pure states |^'„) €Hs(^'Ha, dim Ha = N, such 
that the partial trace Tra|^'n)(^n| — Pn- For pseudopure 
states, the purifications take the form 



l*n> = 



1 - e 
N 



N 



\n) (g) \ipn), 



(35) 



(2) 



where . . . , \^n) is an orthonormal basis of Tia- 7pi 
may be measured interferometrically as a shift in the in- 
terference pattern 



(36) 



by choosing Us ^ e^^iy"-", Vs = ul, Ua = W^™"", 

I 

and Va = (Un) , T being transpose with respect to the 
ancilla basis |(pi),...,|(pjv) and x a variable U(l) phase 

shift m. 

In liquid-state NMR, such a test could in principle be 
implemented as follows. By adding an 'auxiliary' nuclear 
spin qubit pair in the state |0) (g) |0), prepare the total 
input 



IF) = |0) |0) <E) I*, 



(37) 



Apply H(g)H(g)/s(X)/a, H(g)H being Hadamard on each 
auxiliary qubit, followed by the conditional transforma- 
tion 



®\0){0\® Us ®Ua 

|l)(l|®y,®K (38) 
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with Us,Ua,Vs, and Va chosen as above. Finally, apply 
another Hadamard transform on the auxiliary qubits and 
measure the output in the 00 channel. The resulting 
intensity should confirm Eq. 

V. CONCLUSIONS 

We have put forward a kinematic approach to the off- 
diagonal geometric phases of pure states 23!| and nonde- 
generate mixed states 27, 28]. The kinematic approach 
to off-diagonal geometric phases presented in this paper 
emphasizes the path in state space as the primary con- 
cept for the geometric phase and provides a tool to effi- 
ciently calculate off-diagonal geometric phases for arbi- 
trary unitary operators and arbitrary density matrices. 
We have extended the concept of off-diagonal geomet- 
ric phases of nondegenerate mixed states to degenerate 
mixed states. This kind of extention is essential since 
the geometric phase of degenerate mixed states 12j is 
undetermined in the case where the interference visibil- 
ity between the initial and final states vanishes leading to 



an interesting nodal point structure in the experimental 
parameter space that could be monitored in a history- 
dependent manner. On the other hand, the notion of 
off-diagonal geometric phase put forward for nondegen- 
erate mixed state in Refs. |28|j breaks down in the 
case of degenerate mixed states, while the latter case 
does really include the most general quantum states and 
therefore this extention is physically important. Our ex- 
tension makes it possible to obtain information of the ge- 
ometry along the path in state space when the standard 
notion of geometric phase for degenerate density opera- 
tors is undetermined. We have argued that liquid-state 
NMR protocols constitute a possible physical scenario 
for test of off-diagonal geometric phases for degenerate 
mixed quantal states. 
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